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Models 



As long as the commanded values of x produce 0 motion that has a sufficiently 
small value of 0, the approximation given by Eqs. (2.99) or (2.100) is valid and no 
other linearized dynamic relationships are necessary. However, as soon as the 
commanded values of x produce accelerations where the inertial forces (wy and 
the reaction to 19) are a significant fraction of p s - p e , the approximations are 
no longer valid. We must then incorporate these forces into the equations, thus 
obtaining a dynamic relationship between x and 0 that is much more involved than 
the pure integration implied by Eqs. (2.99) or (2.100). Typically, for initial control 
system designs, hydraulic actuators are assumed to obey the simple relationship of 
Eqs. (2.99) or (2.100). 



2.6 Linearization and Scaling 

The differential equations of motion for almost all processes selected for con- 
trol are nonlinear. On the other hand, as will be evident in the next chapter, 
both analysis and control design are far easier for linear than for nonlinear 
models. Linearization is the process of finding a linear model that approxi- 
mates a nonlinear one. Fortunately, as Lyapunov proved over 100 years ago, 
if a small-signal linear model is valid near an equilibrium and is stable, then 
there is a region (which may be small, of course) containing the equilibrium 
within which the nonlinear system is stable. 11 So we can safely make a linear 
model and design a linear control for it such that, at least in the neighbor- 
hood of the equilibrium, our design will be stable. Because a very important 
role of feedback control is to maintain the process variables near equilibrium, 
such small-signal linear models are a frequent starting point for control models. 
Small-signal linearization is discussed in Section 2.6.1. 

An alternative approach to obtain a linear model for use as the basis of 
control system design is to use part of the control effort to cancel the nonlinear 
terms and to design the remainder of the control based on linear theory. This 
approach— linearization by feedback — is popular in the field of robotics, where 
it is called the method of computed torque. It is also a research topic for 
control of aircraft. Section 2.6.2 takes a brief look at this method. Finally, 
some nonlinear functions are such that an inverse nonlinearity can be found to 
be placed in series with it so the combination is linear. This method is often 
used to correct mild nonlinear characteristics of a sensor or actuator that have 
small variations in use. 

The magnitude of the values of the variables in a problem is often very 
different, sometimes so much so that numerical difficulties arise. This was a se- 
rious problem years ago when equations were solved using analog computers, 



In 1949 the Russian scientist Aizerman conjectured that if a certain class of systems were 
stable with any linear gain between two limits, then the nonlinear system with a nonlinear gain 
characteristic that was kept between the same limits would also be stable. Unfortunately, this 
conjecture is not true. 
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and it was routine to scale the variables so that all had similar magnitudes. To- 
day's widespread use of digital computers for solving differential equations has 
largely eliminated the need to scale a problem unless the number of variables is 
very large because computers are now capable of accurately handling numbers 
with wide variations in magnitude. Nevertheless, it is wise to understand the 
principle of scaling for the few cases where extreme variations in magnitude 
exist and scaling is necessary or the computer wordsize is limited. Sections 2.6.3 
and 2.6.4 discuss two kinds of scaling. 

2.6.1 Small-Signal Linearization 

A nonlinear differential equation is one where the derivatives of the state have 
a nonlinear relationship to the state itself and/or the control. In other words, 
the differential equations cannot be written in the form 

x = Fx + Gu 

but must be left in the form 12 

x = f(x,«). 

For small-signal linearization we first determine equilibrium values of x 0 , u 0 , 
that is, values where Xo = 0 = f(x 0 , u Q ). We then expand the nonlinear equation 
in terms of perturbations from these equilibrium values; that is, we let x = 
Xo + 8x and u = u 0 + 8u , so that 

x 0 + 8x = f (x 0 , u 0 ) + ¥8\ + G8u, 

where F and G are the best linear fits to the nonlinear function f(x, u) at x Q 
and u 0 , that is, 

F = I ^ I and G = | ^1 . (2.101) 

Substracting out the equilibrium solution, this reduces to 

8x = F8x + G8u, (2.102) 

which is a linear differential equation approximating the dynamics of the mo- 
tion about the equilibrium point. Normally, the 8 notation is dropped and it is 
understood that x and u refer to the deviation from the equilibrium. 

In developing the models discussed so far in this chapter, we have encoun- 
tered nonlinear equations on several occasions: the pendulum in Example 2.5, 
the hanging crane in Example 2.6, the AC induction motor in Section 2.4, the 



This equation assumes the system is time-invariant. A more general expression would be 
x = f(x, M,r). 
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tank flow in Example 2.20, and the hydraulic actuator in Example 2.21. In 
each case, we assumed either that the motion was small or that motion from 
some operating point was small, so that nonlinear functions were approximated 
by linear functions. The steps followed in those examples essentially involved 
finding F and G in order to linearize the differential equations to the form of 
Eq. (2.102) as illustrated in the following example. 



EXAMPLE 2.22 Linearization of Motion in a Ball Levitator 

Figure 2.34 shows a magnetic bearing used in large turbo machinery. The magnetics 
are energized using feedback control methods so that the axle is always in the center 
and never touches the magnets, thus keeping friction to an almost nonexistent level 
A simplified version of a magnetic bearing that can be built in a laboratory is shown 
in Fig. 2.35, where one electromagnet is used to levitate a metal ball. The physical 
arrangement of the levitator is depicted in Fig. 2.36. The equation of motion of the ball, 
derived from Newton's law, Eq. (2.1), is 



m'x = f m (x, i) -mg, 



(2.103) 



where the force f m (x, i) is caused by the field of the electromagnet. Theoretically, the 
force from an electromagnet falls off with an inverse square relationship to the distance 
from the magnet, but the exact relationship for the laboratory levitator is difficult to 
derive from physical principles because its magnetic field is so complex. However, 
the forces can be measured with a scale. Figure 2.37 shows the experimental curves 
for a ball with a 1-cm diameter and a mass of 8.4 x 10~ 3 kg. At the value for the 
current of 12 = 600 mA and the displacement jci shown in the figure, the magnetic 
force /,„ just cancels the gravity force mg = 82 x 10~ 3 N. (The mass of the ball is 



Figure 2.34 

A magnetic bearing (Photo 
courtesy of Magnetic 
Bearings, Inc.) 
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8.4 x 10~ 3 kg, and the acceleration of gravity is 9.8 m/sec 2 .) Therefore the point (*i, i 2 ) 
represents an equilibrium. Using the data, find the linearized equations of motion about 
the equilibrium point. 

Solution. First we write in expansion form the force in terms of deviations from the 
equilibrium values x\ and i% : 



/ m (x! + 8x, h + 8i) ^ f m (xi, i 2 ) + K x 8x + KiSL 



(2.104) 



The linear gains are found as follows: K x is the slope of the force versus x along the 
curve i = i 2 , as shown in Fig. 2.37, and is found to be about 14 N/m. K t is the change 
of force with current for the value of fixed x — x\ . We find that for i = i x = 700 mA at 
x = x\ , the force is about 122 x 10" 3 N, and at i == 13 = 500 m A at x = x\ , it is about 
42 x 10" 3 N. Thus 



,X?V-: • 

EXAMPLE 



122 x 10- 3 - 42 x 10- 3 80 x 10~ 3 N 



700 - 500 
= 400 x 10" 3 N/A 
= 0.4 N/A. 



200 mA 



Substituting these values into Eq. (2.104) leads to the following linear approximation '^j 
for the force in the neighborhood of equilibrium: - ■:• 

A, = 82 x 10~ 3 + 148x + 0.45/. . 

Substituting this expression into Eq. (2.103) and using the numerical values for mass A || 
and gravity force, we get for the linearized model 



(8.4 x 10" 3 )Jc = 82 x 10~ 3 + 145jc + 0.45* - 82 x 10" 3 . 
Because x = jci 4- 5jc, then jc = 53: . The equation in terms of 8x is thus 
(8.4 x 10-*)8'x = 148x +0.45/, 



(2.105X1 



8x = I6618x + 47.65i, 



which is the desired linearized equation of motion about the equilibrium point. A logical 
state vector is x = [8x 8x] T , which leads to the standard matrices 



'-[1 1] - G =U° 6 ] 



and the control u=8i. 



iations from the 



glAMPLE 2.23 
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Alternate Linearization of the Water Tank 

Repeat the linearization of Example 2.20 using the concepts presented in this section. 
Solution. Equation (2.87) may be written as 

x = f(x,u) (2.106) 

where* 4 h, u 4 Wl ,„and / = -^VpT^* ±w in = ^ -i. _ ^ + J. W/w . 
The linearized equations are of the form 



r approximation 



where 



and 



8x = F8x + G8u, 



2,4/? VPS** - Pa 2Afl VPo - 

[g] = V = _a^ = J_ 



(2.107) 

(2.108) 
(2.109) 

(2.110) 



However, note that some flow is required to maintain the system in equilibrium so that 
Eq. (2.107) is valid; specifically, we see from Eq. (2.87) that 



values for mass 



for h = 0, 



(2.111) 



and the Su in Eq. (2.107) is 8w in , where w in = w ino + 8w in . Therefore, Eq. (2.107) 
becomes 



8h = F8h + G8w in = F8h + Gw in - G-J Po - Pa 



and matches Eq. (2.90) precisely. 



(2.112) 



2.6.2 Linearization by Feedback 

Linearization by feedback is accomplished by substracting the nonlinear terms 
out of the equations of motion and adding them to the control. The result 
is a linear system, provided that the computer implementing the control has 
enough capability to compute the nonlinear terms fast enough. A more detailed 
understanding of the method is best achieved through example. 

To illustrate linearization by feedback, we consider the equation of a simple 
pendulum developed in Example 2.5 [Eq. (2.21)]: 



ml 2 0 + mglsm0 = T c . 



(2.113) 
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If we compute the torque T c to be 

T c = mgl sind +«, 



then the motion is described by 



ml 2 'e = u 



(2.114) 



(2.115) 



Fauation (2 115) is a linear equation no matter how large the angle 6 becomes. 
ESll for purposes of control design because it enables us to use 
Zl aX Schni q ues P £ resulting linear control ^^JS 
of u based on measurements of 0; however, the value of the torque actually 
: „"t to the equipment would derive from Eq. ^^^StSS 
thrP P riaid links this computed-torque approach has led to effective control. 
TTZtlCr^rcnJior the control of aircraft, where the linear models 
change considerably in character with the flight regime. 

2.6.3 Amplitude Scaling 

There are two types of scaling that are sometimes carried out: amplitude scaling 
InSe S Amplitude scaling is usually performed unwittingly by simply 
rickto unhs that make sense for the problem at hand. For the ball levitator, 
Sessing he motion in millimeters and the current in milliamps wotik keep 
he' numbers within a range that is easy to work with. ^^o^siutZ 
developed in the standard SI units of meters, kilograms, and j ampere , but in 
computing the motion of a rocket going into orbit, using kilometers makes 
nZl sens* The equations of motion are usually solved using computer-aided 
dTs nL twaV^ 

sv m it becomes important to scale the problem so that the elements of the 
sSteTec torTave similar numerical variations. A method for accomplishing 
Le 1st sZL for a complex system is first to estimate the maximum value 
for each state element and then Jo scale the system so that each element vanes 

^Tn^erafwe can perform amplitude scaling by defining the scaled vari- 
ables for each state element: If 



then 



x' = S x x, 



= S x x and x' = S x x. 



(2.116) 



(2.117) 



Wethenpick S x to result in the appropriate scale change : su JstituteEq> (2.116) 
Tnd (2.117) into the equations of motion, and recompute the coefficients. 



EXAMPLE 2 .24 Scaling for the Ball Levitator 



Scale the variables for the ball levitator in Example 2.22 to result in units of millimeters 
and milliamps instead of meters and amps. 
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Solution. Referring to Eq. (2.116), we define 

8x' = S x 8x and 8i' = 5,Ji 

such that both S x and 5, have a value of 1000 in order to convert 8x and Si in meters 
and amps to 8x' and 8i' in millimeters and milliamps. Substituting these relations into 
Eq. (2.105) and taking note of Eq. (2.117) yields 

8x f = \6618x' + 47.6— Si". 

Si 

In this case S x = S { , so Eq. (2.105) remains unchanged. Had we scaled the two quantities 
by different amounts, there would have been a change in the last coefficient in the 
equation. 



2.6.4 Time Scaling 

The unit of time when using SI units or English units is seconds. Computer- 
aided design software is usually able to compute results accurately no matter 
how fast or slow the particular problem at hand. However, if a dynamic system 
responds in a few microseconds or if there are characteristic frequencies in the 
system on the order of several MHz, the problem may become ill-conditioned, 
so that the numerical routines produce errors. This can be particularly trou- 
blesome for high-order systems. The same holds true for an extremely slow 
system. It is therefore useful to know how to change the units of time should 
you encounter an ill-conditioned problem. 
We define the new scaled time to be 



r =co Q t (2.118) 

such that, if t is measured in seconds and co Q = 1000, then r will be measured 
in milliseconds. The effect of the time scaling is to change the differentiation 
so that 

dx dx dx 

x== ^7 = ac / x = ^o — , (2.119) 

dt d(x/a> 0 ) dx v ' 

and 

.. d 2 x 2 d 2 x 

X =d^=^d?' (112 °> 

Putting the equation into state-variable form allows a more concise way of 
stating time scaling. For the system described by 

x = Fx + Gw, (2.121) 
we say it is time-scaled, using r = a) 0 t, by the equations 

dx 1 1 

— = — Fx+— Gw. (2.122) 
dx co a co () 
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EXAMPLE 2.25 Time Scaling an Oscillator 



The equation for an oscillator was derived in Example 2.5, For a case with a very fast 
natural frequency w tl = 15,000 rad/sec (about 2 kHz), Eq. (2.23) can be rewritten as 

6 + 15.000 2 0 = 10 6 • T c . 

Determine the time-scaled equation so that the unit of time is milliseconds. 

Solution. The value of co 0 in Eq. (2.118) is 1000. Equation (2.80) shows that 



and the time-scaled equation becomes 



— + 15 2 -0 = T c 
dr 2 



In practice, we would then solve the equation 

9 + 15 2 9 = T c (2.123) 



and label the plots in milliseconds instead of seconds. 

In state-variable form with a state vector x = [0 6] T , the unsealed matrices are 



Applying Eq. (2.122) results in 



which yields state-variable equations which are equivalent to the scaled system in 
Eq. (2.123). 
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SUMMARY 

Mathematical modeling of the system to be controlled is the first step in 
analyzing and designing the required system controls. In this chapter we 
developed models for representative systems. Important equations for 
each category of system are summarized in Table 2.1. 

An alternative way of expressing the differential equations that character- 
ize the model of a linear system is the state-variable form, 

x = Fx + Gw, 
y = Hx + Ju. 



TABLE 2.1 Key Equations for Dynamic Models 



important Laws Associated 

System or Relationships Equations Equation Number 

Mechanical Translation motion F = ma (2.1) 
(Newton's law) 

Rotational motion M = la (2.14) 

Motion of nonrigid bodies mx + bx + kx = F (2.10) 

State-variable form Linear system x = Fx + Gw (2.39) 

y = Hx + /w (2.40) 

Electrical Operational amplifier (2.47), (2.48) 

Electromechanical Law of motors F = Bli (2.54) 

Law of the generator e(t) = Blv (2.57) 

Torque developed in a rotor T = K t i a (2.61) 

Back emf Voltage generated as a result e = K e 6 m (2.62) 
of rotation of a rotor 

Heat flow Heat-energy flow q=y(7i- 72) (2.73) 

Temperature as a function t = (2.74) 

of heat-energy flow 

Specific heat C = mc v (2.75) 

Fluid flow Continuity relation m = w m - w out (2.81) 
(conservation of matter) 

Force of a fluid acting / = pA (2.83) 
on a piston 

Effect of resistance w = -j^ (pi - p 2 ) Xfot (2.84) 
to fluid flow 
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Equations in state-variable form are conducive to solution by computer 
packages that were developed especially for matrix equations (e.g., Matlab). 
In Section 2.2 we briefly introduced the state-variable form; it will be ex- 
plored in more depth in Chapter 7. 
. Linearization and scaling (Section 2.6) are methods by which certain com- 
plications of dealing with differential equations can be minimized. In 
linearization, nonlinear differential equations are approximated by linear 
ones by either (1) considering a small-signal linear model that is acurate 
near an equilibrium, or (2) linearization by feedback, or (3) introducing an 
inverse nonlinearity. Scaling of variables results in numerical values that 
fall within a narrow-enough range of magnitude to minimize errors and 
allow for ease of computation. 

Review Questions 

1. What is a "free-body" diagram? 

2. What are the two forms of Newton's law? 

3. Why is it convenient to write equations of motion in the state-variable form? 

4. For a structural process to be controlled such as a robot arm, what is the meaning of 
"collocated control"? "Noncollocated control"? 

5. When, why, and by whom was the device named an "operational amplifier"? 

6. What is the major benefit of having zero input current to an operational amplifier? 

7. State Kirchoff s Current Law. 

8. State Kirchoff's Voltage Law. 

9. Why is it important to have a small value for the armature resistance, K„, of an 
electric motor? 

10. What are the definition and units of the electric constant of a motor? 

11. What are the definition and units of the torque constant of an electric motor? 

12. Give the relationships for (a) heat flow across a substance and (b) heat storage in a 
substance. 

13. Name and give the equations for the three relationships governing fluid flow. 

14. Why do we approximate a physical model of the plant (which is always nonlinear) 
with a linear model? 

Problems 

Problems for Section 2.1 

2.1. Write the differential equations for the mechanical systems shown in Fig. 2.38. 

2 2 Write the equations of motion of a pendulum consisting of a thin, 2-kg stick of 
length / suspended from a pivot. How long should the rod be in order for the 
period to be exactly 2 sec? (The inertia / of a thin stick about an endpoint is 
\ml 2 . Assume 6 is small enough that sin0 = 9.) 
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